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Influence of steps on the tilting and adsorption dynamics of ordered pentacene films 

on vicinal Ag(lll) surfaces 
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Here we present a structural study of pentacene (Pn) thin films on vicinal Ag(lll) surfaces by 
He atom diffraction measurements and density functional theory (DFT) calculations supplemented 
with van der Waals (vdW) interactions. Our He atom diffraction results suggest initial adsorption 
at the step edges evidenced by initial slow specular reflection intensity decay rate as a function of 
Pn deposition time. In parallel with the experimental findings, our DFT+vdW calculations predict 
the step edges as the most stable adsorption site on the surface. An isolated molecule adsorbs as 
tilted on the step edge with a binding energy of 1.4 eV. In addition, a complete monolayer (ML) 
with pentacenes flat on the terraces and tilted only at the step edges is found to be more stable than 
one with all lying flat or tilted molecules, which in turn influences multilayers. Hence our results 
suggest that step edges can trap Pn molecules and act as nucleation sites for the growth of ordered 
thin films with a crystal structure similar to that of bulk Pn. 

PACS numbers: 82.45.Mp, 68.55.-a, 73.21.Ac, 71.15.Mb, 68.47.Fg 



Pentacene thin films are still being studied heavily due 
to their potential electronic device applications and being 
a model system for organic semiconductor film studies, [lj- 
3] On many different metal surfaces, Pn is reported to 
form lying down films with similar unit cell structures. [4]- 
[bjj In case of the second or higher layers however, even 
on the same metal surface, different growth mechanisms 
and/or dynamics have been reported. [J, |8|, |9|, ll5H17| On 
Ag(lll), while some reports suggest a bilayer film for- 
mation, where an ordered second layer with the symme- 
try of the Ag(lll) surface forms on top of a disordered 
(2D gas phase) first layer at room temperature, some 
others suggest formation of bulk like Pn structures im- 
mediately after the first ML).[4j, |l6|, ll7| In addition, in 
our earlier works we had shown that an ordered Pn mul- 
tilayer could only be formed on a relatively high step 
density Ag(lll) surface and the multilayer order can be 
improved by increasing the kinetic energy (KE) of the 
molecules, by seeded supersonic molecular beam (SSMB) 
deposition. [131] This observation was attributed to step 
edges acting as nucleation centers for tilted molecules 
in the multilayers resulting in a step flow growth mecha- 
nism. Such growth dynamics was also observed for differ- 
ent planar organic molecules on vicinal copper and gold 
surfaces by different groups. [18H20} In addition, high en- 
ergy of the molecules was suggested to result in a local 
annealing effect which in turn helps formation of ordered 
multilayers at low substrate temperatures. These exper- 
imental results were in line with our recent theoretical 
work on Pn on flat Ag(lll) surface where a multilayer 
structure with tilted Pn molecules was found to be more 



stable than one with flat lying molecules. |21( However 
during the above mentioned experimental studies surface 



step density was changed arbitrarily by using a miscut 
Ag(lll) surface and the molecular KEs were not mea- 
sured by a time of flight (TOF) method but only esti- 
mated by using a theoretical formula. 13| The necessity 
to test the stable configuration and the contradictions in 
the literature have stimulated the present experimental 
work on Pn film growth on Ag(ll,12,12) surface comple- 
mented by density functional theory (DFT) calculations 
to explore the effect of surface step density and molecular 
energy on the Pn film growth with the aim to improve 
our understanding of such systems. 

Pn films were grown by using SSMB deposition and 
characterized with a low energy atom diffraction (LEAD) 
apparatus both of which were detailed before. [12u22| KE 
of the molecules was controlled by changing the noz- 
zle temperature of the SSMB source and the carrier gas 
(CG). A commercial quartz crystal microbalance (QCM) 
and a quadruple mass spectrometer were implemented 
on this apparatus for allowing molecular flux and TOF 
measurements. In addition, a gold covered electrode of a 
custom design QCM can be used as the active substrate 
to simultaneously measure He specular intensity decay 
and record resonance frequency shift of the substrate to 
estimate film coverage during deposition. [23j Pentacene 
KEs, determined by TOF measurements, for standard 
SSMB source conditions were 2.2 eV for He CG and 0.2 
eV for Kr CG. For LEAD measurements a monoenergetic 
He beam, with energy of 14 meV (incident wave vector 
fcj=5.13 A" 1 ) and velocity dispersion of 2%, was scat- 
tered for the sample surface and the angular distribution 
of the scattering intensity was recorded. For reciprocal 
space mapping the diffraction spectra, recorded at 80 K 
surface temperature along different crystallographic di- 



rections, were converted to momentum by using the equa- 
tion AA'n = ki(sin6f — sin#j) where Aifii is the parallel 
momentum transfer, and 9i and 9f are the incident an- 
gle and the detector angular position, respectively. The 
reference "flat" Ag(lll) and the Ag(ll,12,12) substrates 
were purchased from Mateck GmbH with an orientation 
accuracy of <0.1° and had specular reflection intensity 
values of 60 % and 45 % of the incident beam respec- 
tively Ag(ll,12,12) has a miscut of 2.31° with respect 
to Ag(lll) surface normal, with (111) terraces of 24 rows 
(58 A) wide with step edges parallel to [Oil] direction and 
perpendicular to [211] direction. 
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FIG. 1. Specular reflection intensity decay curves as func- 
tion of Pn coverage, (a-d) Insets show the initial decay region. 
(ki—5.13 A -1 , incidence angle is 56° for Ag(lll) and 59° for 
Ag(ll, 12,12), deposition temperature 200 K, on Ag(ll, 12,12) 
molecular beam direction was parallel to step edges) In (e) 
specular reflection intensity decay curve and resonance fre- 
quency shift from the gold electrode surface of a quartz crystal 
during Pn film deposition is shown. 

In FigfT] we show He specular reflection intensity de- 
cay curves as a function of film coverage, at 200 K sub- 
strate temperature, for different surfaces and deposition 
energies. The coverage values were deduced from the 
specular reflection recovery peaks, with the assumption 
that these peaks correspond to ML completion. In fact, 
a similar recovery peak was measured on Au(lll) elec- 
trode surface of a quartz crystal at a resonance frequency 
shift (2.9 Hz) corresponding to a deposited ML of Pn as 
shown in FigJlJ. When the initial part of these curves 
are examined it can be seen that, regardless of the type 
of the CG, on Ag(ll,12,12) there is a slow decay up to 
about 0.12 ML that is due to initial Pn adsorption to 
step edges. In fact, the He reflectivity is very sensitive to 
molecular adsorbed species [24j and this slow decay can 
be attributed to the overlapping of Pn cross sections (a) 
with that of silver step edges resulting in a decrease of the 




FIG. 2. Helium diffraction scans of mono (a) and multilayers 
(b) of Pn on Ag(ll,12,12) with different KEs. Comparison 
of multilayer diffraction scans on Ag(ll,12,12) and Ag(lll) 
grown with high KE are shown in (c). (ki—5.13 A -1 , incidence 
angle is 56° for Ag(lll) and 59° for Ag(ll,12,12)) 



effective cross section (cr e ff ) of Pn molecules. However on 
Ag(lll) due to lack of steps, Pn molecules mainly adsorb 
on the terraces and present a larger a e s to incoming He 
atoms which results in a faster decay of specular inten- 
sity. On the 24 rows wide terraces of Ag(ll,12,12), 8 Pn 
rows with a 6x3 unit cell are expected to adsorb at full 
coverage [lj 



2l| and a single Pn row is necessary to sat- 



urate the steps, corresponding to a coverage 1/8 (0.125) 
ML. For higher coverage, Pn starts to fill terraces and 
the decay is faster explaining the observed specular in- 
tensity trend. In addition, when deposition is performed 
at 80 K substrate temperature, the slow decay region 
is completely missing due to reduced diffusion length of 
Pn molecules which prevent them to reach to the step 
edges. Regardless of the surface step density, the reflec- 
tivity decay curves for coverages of >0.3 ML show that 
deposition with Kr yields a smoother ML as indicated 
by the higher specular reflection intensity level. This can 
be due to disordering of the Ag step edges and terraces 
caused by the high energy of the Pn molecules deposited 
with He CG. In fact for Ag(lll) surface, energies needed 
to detach an atom from a step edge (0.49-0.76 eV) or to 
diffuse an atom along the step (0.29-0.34 eV)[25|, |26[ are 
all much lower than Pn KEs deposited with He. Whereas 
in case of Kr CG, Pn molecules do not have enough en- 
ergy to activate these roughening processes. In fact Pn 
cr c ffS on Ag(ll,12,12) (Fig [TH), calculated by maximal 
attraction model proposed by Comsa 24| , increases from 
(7 K r = 102 ± 3 A 2 to ctho = 156 ± 4 A 2 for deposition 
with Kr and He respectively. When deposition was per- 
formed with a He beam perpendicular to step edges, cthc 
further increases to 328 ± 2 A 2 , which indicates that step 
edge roughening is more effective in case of pcrpendicu- 



These values can be compared with the 
124A 2 [2l-119 A 2 [li. Af- 



lar deposition 

Pn vdW cross section, <7 v dw 

ter the completion of the first ML, specularity decreases 

much faster in the case of Kr CG on both surfaces. This 

suggests that though the ML is more ordered/smooth 

with Kr, an ordered second layer does not grow on it. 

Diffraction scans, obtained from the MLs and the mul- 
tilayers grown with different KEs and on different sub- 
strates, shown in Fig. [2] confirm the above mentioned 
mechanism: On Ag(ll,12,12) while the ML diffraction 
peaks and specularity are more intense with Kr CG, no 
multilayer diffraction peaks could be observed. A peak 
width analysis indicates that average ML domain size 
with He beam is about 50 A whereas with Kr it is about 
70 A which means that with both CGs the domains ex- 
tend as wide as the terrace width of the Ag(ll,12,12) 
surface. With He, however, multilayer diffraction peaks 
which are consistent with previously reported[12j, ll3( Pn 
multilayer structure could be observed. Finally when the 
multilayer structure on Ag(lll) and Ag(ll,12,12) are 
compared it can be seen that, though less intense, the 
diffraction peaks are slightly narrower on Ag(ll,12,12), 
corresponding to domain size of 17 A on Ag(lll) and 24 
A on Ag(ll,12,12). The reason for low specular reflec- 
tion and diffraction peak intensities on the vicinal surface 
may be the higher level of diffuse scattering due to step 
edges, as discussed above. 

To gain further insights on the system, we performed 
DFT calculations. Pure DFT methods underestimate 
the adsorption energy of Pn molecules on Ag(lll). 
Among these, GGA-PW92 functional [H] gives still low 
but slightly better binding energies and bond lengths for 
isolated and ML cases on Ag(lfl) surface. |21| DFT cal- 
culations were performed using VAS PI30I with projec- 
tor augmented waves (PAW) method [3l|, |32J. We em- 
ployed Grimme's semiempirical forcefield approach[33| to 
include long range vdW corrections to better describe the 
interaction of Pn molecules with Ag substrate. We used 



TABLE I. DFT+vdW results of Pn on vicinal Ag(lll) surface 
for different coverage models. E t is the relative supercell total 
energy. Pn-step edge bond length is denoted by I. Labels 
denote the layer numbers for interlayer separations, di, di-2 
and for tilting angles ai, Q2, respectively. 





Model 


Et(eV) 


1(A) 


rfi(A) 


di- 2 (A) 


ai(°) 


a a (°) 


Single 


flat@bridge 
C@step 


0.00 
-0.47 


2.55 


2.92 


- 


0.0 
25.4 


- 


1ML 


flat 
flat2 a 


0.00 
-0.51 


2.63 
2.49 


2.94 
2.94 


— 


0.0 

0.0 


— 


2ML 


flat-tilt 
flat2 a -tilt 


0.00 
-0.49 


2.54 
2.50 


2.77 

2.77 


2.56 
2.56 


0.0 
0.0 


17.0 
-13.1 
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molecules flat on the terrace, tilted at the step edge. 



b) 



FIG. 3. (a) Potential energy profile of single Pn on Ag(233) 
vicinal surface along [211] direction (A,B,C configurations 
shown in the corresponding insets), (b) the minimum energy 
configuration of an isolated Pn at the step edge on Ag(455) 
surface (the charge density is presented in the inset). 



a 4 layer slab to model the vicinal Ag(lll) surface with 
two different terrace sizes. The width of the terraces on 
Ag(233) and Ag(455) supercells measure 10.0 A and 20.0 
A spanning six and ten surface Ag atoms along [211], 
respectively. Ag(233) supercell has been used to study 
the potential energy profile of an isolated Pn as schemat- 
ically shown in FiglHK- Ag(455) slab for one and two ML 
coverages allows to place three Pn molecules one being 
at the step edge as depicted in FiglU We have relaxed 
all probable configurations based on minimization of the 
forces requiring them to be smaller than 0.1 eVA -1 on 
every atom. In order to elucidate the effect of van der 
Waals contribution on Pn-metal interaction that charac- 
terizes the minimum energy structures, we carried out 
these calculations both with and without dispersive en- 
ergy corrections. 

Pn molecules give noticeably stronger binding at the 
step edges relative to that on the flat terraces as seen 
in FigUJi which shows the diffusion barrier profile along 
[211] obtained with DFT. Pn at the step edge bends along 
its long axis such that the middle part of the molecule 
gets slightly more closer to the surface from where we 
measure Pn-step edge bond length (column I in Table [TJ) . 
The curvature is similar to the one observed for Pn on 
Cu(lll) surface. [34[ DFT+vdW reproduces this curva- 
ture that indicates a strong adsorption. Binding energies 
at the step edge (case C in Fig(3K) and at the bridge site 
(molecular center on a Ag-Ag bond) on the terrace are 
1.40 and 0.92 eV, respectively. These were much lower 
with pure DFT. For instance, DFT predicts the binding 
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FIG. 4. Relaxed geomtries of 2ML Pn on Ag(455) for (a) 
flat-tilt and (b) Hat2-tilt cases. Densities of states (DOS) 
are presented for the flat2 and flat2-tilt structures where the 
shaded regions are Pn contributions. 



energy at the step edge to be 0.62 eV with no bending. 

For the first ML, we have considered tilted, flat and 
flat2 models as initial structures for geometry optimiza- 
tion. In tilted case all molecules are tilted around their 
major axes. Flat model is as the first layer of FigHk 
where all Pn molecules lie parallel to the Ag surface 
plane. Flat2 case differs from the flat one in only that the 
Pn at the step edge is tilted as in the most stable case C of 
an isolated molecule (see the first layer of FigHb). Pure 
DFT predicts tilted formation to be the minimum energy 
configuration while DFT+vdW relaxes tilted case to flat2 
one. Moreover, the measured a e ss provide a further sup- 
port to flat2 model because Pn covers completely a step 
edge and the tilting reduce the a to cr v dw cos(25.4) = 112 
A which is at most 10% greater than <7Kr- In the flat 
model, the cr v dw is increased by the step edge cross sec- 
tion, so, apparently the chc supports the flat case. But 
the experimental evidence that He deposition causes an 
increased disorder and further calculations that exclude 
the presence of other configurations with minimal energy 
similar to flat 2, leaves flat2 as the only possibility. 

Then we considered a second ML which can be tilted or 
flat as an initial configuration. For 2ML Pn on Ag(455), 
pure DFT predicts tilted— tilted model to be the lowest 
energy structure. When the geometry optimization re- 
peated with dispersive corrections, tilted molecules on 
the terrace become flat leading to flat2-tilt configura- 
tion (FigHjD). Therefore, second ML is tilted while first 
layer lies flat similar to that of the 1ML case. This is in 
good agreement with the proposed model for bilayer Pn 
on Cu(lll) from STM/STS measurements. Q This in- 
dicates the role of the step edge in bulk-like formation of 
Pn molecules beyond the second ML. As a trend, energy 
differences between flat and flat2 cases for both 1ML and 
2ML coverages follow from that of single molecule adsorp- 
tion cases, i.e. flat@terrace and C@step as presented in 
Table Q] Moreover, a comparison of DOS structures for 
Ag(455) in Fig|4]with those obtained for Ag(lll) without 
steps [21| shows that the difference between the HOMO 
level of Pn and the Fermi energy of the surface decreases. 
This suggests an enhanced charge carrier injection at Pn- 
metal interface in favor of the vicinal surface. 

In summary, by using He diffraction and QCM tech- 
niques simultaneously and performing DFT+vdW calcu- 



lations we have, unambiguously, shown that step edges on 
Ag(lll) surfaces act as nucleation centers for formation 
of bulk-like Pn thin films on MLs composed of fiat lying 
molecules. This mechanism can be considered to be a 
rather general one since our results are in line with those 
reported on different metal surfaces such as Au and Cu. 
Our results highlight the importance of a) vdW forces for 
systems with relatively weak interactions like Pn/metal 
interfaces and b) the effects of metal substrate surface 
morphology on the crystal and electronic structure of or- 
ganic semiconductor films which have very important im- 
plications on the performance of electronic devices that 
employ these metal/organic interfaces. 
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